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Abstract 

Wilf posed the following problem: determine asymptotically as n ^ oo the prob- 
ability that a randomly chosen part size in a randomly chosen composition of n has 
multiplicity m. One solution of this problem has been given by two of the authors Q. 
In this paper, we study this question using the techniques of generating functions and 
singularity analysis. 
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1 Introduction 



Let n be a positive integer. A composition of n with p parts is a solution of the equation 
= Ki + /t2 + ■ ■ ■ + in positive integers ki, K2, . . . , Hp. We shall write k = K2, ■ ■ ■ , i^p) 
to symbolize the composition. For example, there are 16 compositions of 5, namely 

(5) (4,1) (1,4) (3,2) 

(2,3) (3,1,1) (1,3,1) (1,1,3) 

(2,2,1) (2,1,2) (1,2,2) (2,1,1,1) 

(1,2,1,1) (1,1,2,1) (1,1,1,2) (1,1,1,1,1). 

The terms called the parts of the composition. The multiplicity of a part 

size is the number of parts with that size. For example, in the composition (1, 1, 1,2) the 
multiplicity of 1 is 3 and the multiplicity of 2 is 1. A partition of n with p parts is a solution of 
n = Xi + \2 + - ■ ■ + Xp with Ai > A2 > ■ ■ • > Ap. In [|l| Corteel, Pittel, Savage, and Wilf proved 
that for every fixed m > 1, the probability that a randomly chosen part size in a random 
partition of n approaches l/(m(m + 1)) as n — > 00. Wilf then posed the corresponding 
problem for compositions: determine asymptotically (as n 00) the probability that a 
randomly chosen part size in a randomly chosen composition of n has multiplicity m. One 
solution of this problem has been given by two of the authors 0]. In this note, we address 
the same question using generating functions and singularity analysis. 

It is well known that there are 2"~^ compositions of n. One way to arrive at this result 
uses generating functions. The generating function for compositions with p parts is 

(z + z'^ + z^ + ---y-- 



1 - 

and summing over p we have the generating function for all compositions: 



G{z) 



00 



^\l-zj l-2z 
p=i 



The coefficient of in the expansion of G{z), denoted by [z"]G'(z), is the number of compo- 
sitions of n, and clearly [2;"]G(2;) = 2"~^. From an analytic point of view, the fact that the 
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number of compositions of n is asymptotically (as well as exactly) 2" ^ is associated with 
the fact that the generating function is a rational function for which the pole nearest the 



problem uses the same principle. 

We shall use the following notation. The probability of the event A is denoted by P(v4), 
and the expected value of a random variable X is denoted by E(X). The natural logarithm 
and base 2 logarithm are denoted by log n and logg n, respectively. 

To state the problem more precisely, suppose that a composition k is selected uniformly 
at random from the set of all 2"~^ compositions of n. Then out of the set of part sizes in 
K, a part size k is chosen uniformly at random. Let aI"^^ denote the event in which k has 
multiplicity m. For example, inspection of the 16 partitions of 5 shown above yields 



and otherwise f[Ai^ ) = 0. The object is to determine F{An ) asymptotically as n — oo. 
We shall find that F{A^^) tends to at the rate 1/ logra. It then turns out that logn-F{A^^) 
does not have a hmit, but oscillates about the value 1/m as n oo. 



The answer to Wilf's question is given in the following theorem, first proved in P|. 

Theorem 1. Let An™^ be the event in which a randomly selected part size in a randomly 
selected composition of n has multiplicity m. Then 



origin (the only pole in this case) is simple and located at z 



|. Our solution of Wilf's 




2 Results 



logn ■ P(A(r)) = (1 + o(l)) - + F({log2 n}) 




n 



oo 



where {a} 



= a — 



[a\ is the fractional part of a and 




with r denoting the gamma function. 
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Using well-known facts about the gamma function (r(l-|-2;) = zr{z) and r{z)r{l — z) — 
tt/ sin(7r2;)), we obtain 

= —,J2{ • u ^ cos(27rpx - 0^), 
ml \smn(pQ;) / 

where a — 27r^/log2 and (f)p is the argument of r(l + i27rp/log2). This series converges 
quite rapidly, and its sum may be approximated by the first term. But even the first term 
is quite small since 2(a/ sinha)^/^ ^ 10~^. Thus for large n one finds that logn • ¥(A^^) is 
quite close to 1/m, but there is a residual dependence on {log2n}. In the treatment given 
here using generating functions and singularity analysis, the proof of Theorem 1 will reduce 
to a well-known calculation after we have established the appropriate sequence of lemmas. 

Let K be a composition of n. Then 'P(k) will denote the set of distinct part sizes in k, 
and J\4m{K) will denote the set of part sizes of k that have multiplicity m. 

Lemma 1. For a random composition of n, 

Proof. Let Gk{z,w) be the two-variable generating function in which [z"w'^]Gk{z,w) is the 
number of compositions of n in which k has multiplicity m. To construct such a generating 
function, we first note that the contribution made by compositions with p (not necessarily 
distinct) parts is 

(z + z^ + --- + WZ'' + z''+^ + i^— + {w- ^)z^ 

\1 — z 
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Thus 



p=i ^ ' 

z + (w - '\)z^{\ - z) 



l-2z-{w-l)z''{l-z) 
1-z 

l-2z- {w-l)zf'{l- z) 
1-z 



- 1. 



l-2z + z^{l -z)- wz^{l - z) 
Since there are 2"~^ compositions of n, we then have 

fik e MM) = T^iz^'w'^]- — 177^—4 TT. r 

V my)) 2n-lL ^l_2z + Z^{1- Z) -WZ^{1- Z) 

_ 1 ^fem(^ _^)m+l 

~ 2^^^ \l-2z + z\l-z))^+^' 
as claimed. □ 

Lemma 2. The polynomial l — 2z + z^{l — z) has precisely one zero z = p satisfying \z\ < 1. 
This zero is given by 

1 1 ^ ^ ^ ^ 



2 2^+2 V22'= 

For allk > 1, 



(-^) < (2^ < (-2^) ■ 



(2p) 

Proof. For the first part, simply observe that if 2; = e*^ then |1 — 2z\'^ = 5 — 4cos^^ and 
|1 — 2;p = 2 — 2 cos so |1 — 22;| > |1 — 2;| for all z with |2;| = 1. Apply Rouche's theorem. To 
get the approximate location of p, write p = | + e and substitute into 1 — 2p + — p) = 0. 
This yields 

1 k-1 k(k-3) o 

€ — 1 e H € + ■ ■ ■ 

2k+2 ^ 2*^+1 2^= ' 
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and thus the stated result by iteration. Next we prove 

1 1 

< P< 7: + 



A simple calculation shows that Q{x) — 1 — 2x + x''{l — x) decreases on (0, 1). Set a — 
1/(2 - 2-(*^+^)) and 6 = I + 2-(*^+^. Then we find that Q{a) > and Q{b) < 0, so a < p < 6. 
Then since (1 + x)" < exp{nx) for x > —1, we have 



(2^<(l-2^) <^-p(-2^) (2^^ (l + 2-^)n >^"P(-^) 



□ 



We shall show that the number of distinct part sizes of a random composition of n 

satisfies ~ log2n with probability 1 — o(l). The underlying probabihstic considerations 
are given in the following lemma. 

Lemma 3. Let X — Y^Ij where (Ij) are indicator random variables. Suppose that P(/fc) = 
Pk,n = Pk- If ci and b are chosen so that both J2j<ai^ ~ Pj) ^'^^ Ylj>bPj '^''"^ then 

F{a<X <b)> 1 -o(l). 

Proof. For all a < 6 we have 

F(a<X <b)^l- F{{X <a}U{X >b})>l- ¥{X < a) - F(X > b). 
Now, denoting for simplicity a set and its indicator by the same symbol, 

nx > 6) < p ( u /, ) < j]p(/,) = J2pv 

\j>b ) j>b j>b 

P(X < a) < P ( U A < 5^(1 - P(/,)) = ^(1 - p,). 

\j<a / j<a j<a 

Hence, if a and b are chosen so that both X^j<(j(l — Pj) and J2j>bPj ^i^) 

¥(a<X <b)> 1 -o(l). 

□ 



and 
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Lemma 4. Let k be a random composition of n. As n oo the number of distinct part 
sizes satisfies ~ log2 with probability 1 — o(l). 

Proof. As a special case of Lemma 1, the probability that k has multiplicity in the random 
composition n is 



Hence 



From Lemma 2, the rational function (1 — 2;)/(l — 2z + z^{l — z)) is analytic for \z\ < 1 
except for a simple pole at z = p ^ The residue is —(1 — p)/{2 + {k + l)p'' — kp'^). By 
standard arguments H, §5.2], 



By Lemma 2, 



^ h-2z + z^{l-z) V2 + (A; + l)p^-A;pfc-V P"+' 



1^ + 0^* 



2 + (A; + l)p'= - fcp'^-i 4 V 
Hence we have 



2„-iL Ji_2;2 + ^fc(i_;2) (2p)"+i V V2'' 
Using the general bound from Lemma 2 

f n\ 1 ( n 

exp -7^ < 777^ < exp 



2V (2p)" ^ ^ ^ V 2'=+2 

we see that 

e = 1 - , \ (l + o(X 

+ 1 r n + 1] ^ f k 
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so that letting b — [log2 n\ + log log n we get 



k>b 

Similarly, 



F{keV{K))>l-exp\-^\(l + 0^'' 



2^+1 I \ \ 2*^ 



that is 



1 - P(*; e BM) < exp |-^i±i} (l + O (A)) . 



Consequently, for any positive a, 



L<fe<a l<fc<a ^ 

0<r<a ^ ■' 

<C5]exp{-(r + l)^} 



= C 



r>0 

exp{-(n + l)/2"+2} 
1 -exp{-(n+ l)/2"+2}^ 



and thus 



l<fe<a 



provided a < [log2nJ — loglogn. Hence, by Lemma 3 applied to Ik = {k E T>{k)}, \D{k)\ ~ 
log2n with probability 1 — o(l). □ 



Given a random composition k, the probability that a randomly selected part thereof has 
multiplicity m is Lemma 4 greatly simplifies the basic problem. Since so 

doing amounts to the neglect of a set of compositions with total probability measure o(l), we 
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may assume that as n — > oo the randomly selected composition k satisfies ~ loggn. 

Thus P(4"^) ~ E(|M^|)/log2n. 

Now we wish to study the asymptotic behavior of P(A; e /Arn{i^)), with the aim of 
estimating 

HlMml) = Y.^{k e MM). 

k 

Lemma 5. The expected value of \A4rn\ is given by 

E{\M^\) = (1 + 0(1))^ V2-^™exp(-n/2'=). 

ml ^-^ 

k 

Proof. As we found in Lemma 1, the relevant generating function is 

^(^) = ^gSfrli)' ^here P{z) = z'-{l - zr^\ Q{z) = 1 - 2z + z'{l - z). 

Recall that Q has a simple zero at z = p ~ ^ and no other zeros in {z : |-2|<1}. In a 
deleted neighborhood of p, we have the Laurent expansion 

„/ X m+1 ~" 



s=0 



The asymptotic behavior of [z"']P{z)/Q"^~^^{z) is governed by the principal part, more specif- 
ically by the r = m + 1 term. In view of 



[^"](1 -^)-" = 

a simple calculation gives 



n + a — 1 
n 



' \ m J {-pQ'{p))m+i {2p)n \ \n 

Set 

logn — log log n — log(4(m + 1)) 



q{n) = 



log 2 
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and note that if A; < q{ri) then 2*^ < n/{A{m + 1) logn), so 

< n"*exp ^2fc ) ^'"exp(— (m + 1) logn) = — . 



(2p) 
Hence we have 

nk & MM) ^ o {^''^'' 

k<q{n) 



n 



n — > oo. 



In view of the fact just noted, in estimating ^i^F{k G A4m), we can hmit ourselves to 
cases where k > q{n). In that case 

jk^-M-^){^-o{^)). «... 

Now 

p(^) ^ 2-^-2-(-+^) (l + O (i^)) and pg'(p) = -1 + O (i^) , 

so 

P(* e A^^-c)) = (l + O ) ^ 2-('«)» exp (-^) . 

It is now evident that the contribution to the sum ^^IP(/c G Aimin)) from those terms 
with k > log2n + loglogn is o(l), so there are O(loglogn) terms in the sum that make a 
nontrivial contribution. Thus the bound on the error for an individual term suffices to give 
the correct asymptotic result for the sum. Replacing /c + 1 by /c in the sum, we have the 
stated result. □ 

Proof of Theorem 1. The computational problem that remains is the asymptotic evaluation 
of 



k=l 



^2-'^™exp(-n/2'=). 



Problems of this kind occur frequently in probability theory and the analysis of algorithms, 
and now there are different methods available for their study, and these methods are described 



10 



in several excellent references 0, chapter 7]. We sketch an approach due to N. G. de Bruijn, 
which is described in pp. 131-134] and elsewhere. A special case (m = 1) of the above 
sum is treated in 0. The starting point is Mellin's formula 

exp(— to) = — : / w~^T{z)dz, w, 0" > 0. 

Substituting this representation (with cr = m — |) for exp(— ri/2''') and using uniform con- 
vergence, one obtains 



y2-'^™exp(-n/2'=) = — T '~ ^i-HM rf^. 
Then by the residue theorem, 

^j;2-'=-exp(-n/2'=) 



ml 

k=l 



J (m-l)! + 2Re e-^^^P^"^^" T (m + 



^1^) Hi + o(i)). 



m!log2 I \ log 2 

The stated result for logn ■ P(yli™''') follows. □ 
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